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Abstract 

We study propagation algorithms for the conjunction 
of two AllDifferent constraints. Solutions of an 
AllDifferent constraint can be seen as perfect matchings 
on the variable/value bipartite graph. Therefore, we inves- 
tigate the problem of finding simultaneous bipartite match- 
ings. We present an extension of the famous Hall theorem 
which characterizes when simultaneous bipartite matchings 
exists. Unfortunately, finding such matchings is NP-hard in 
general. However, we prove a surprising result that finding a 
simultaneous matching on a convex bipartite graph takes just 
polynomial time. Based on this theoretical result, we pro- 
vide the first polynomial time bound consistency algorithm 
for the conjunction of two AllDifferent constraints. We 
identify a pathological problem on which this propagator is 
exponentially faster compared to existing propagators. Our 
experiments show that this new propagator can offer signifi- 
cant benefits over existing methods. 



Introduction 

Global constraints are a critical factor in the success of con- 
straint programming. They capture patterns that often oc- 
cur in practice (e.g. "these courses must occur at different 
times"). In addition, fast propagation algorithms are associ- 
ated with each global constraint to reason about potential so- 
lutions (e.g. "these 4 courses have only 3 time slots between 
them so, by a pigeonhole argument, the problem is infeasi- 
ble"). One of the oldest and most useful global constraints is 
the AllDifferent constraint ( [Lauriere 1978| l. This speci- 
fies that a set of variables takes all different values. Many 
different algorithms have been proposed for propagating 



the AllDi fferent constraint ( |Regin 1994[|Leconte 1996 



Puget 1998j ). Such propagators can have a significant impact 
on our ability to solve problems ( [Stergiou & Walsh 1999| l. 

Problems often contain multiple AllDifferent con- 
straints (e.g. "The CS courses must occur at different times, 
as must the IT courses. In addition, CS and IT have several 
courses in common"). Currently, constraint solvers ignore 
information about the overlap between multiple constraints 
(except for the limited communication provided by the do- 
mains of common variables). Here, we show the benefits of 
reasoning about such overlap. This is a challenging problem 
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as finding a solution to just two AllDifferent constraints 
is NP-hard ( |Kutz etal. 2008 ) and existing approaches to deal 
with such overlaps require exponential space (Lardeux et al. 
2008| l. Our approach is to focus on domains that are ordered, 
as often occurs in practice. For example in our time-tabling 
problem, values might represent times (which are naturally 
ordered). In such cases, domains can be compactly rep- 
resented by intervals. Propagation algorithms can narrow 
such intervals using the notion of bound consistency. Our 
main result is to prove we can enforce bound consistency on 
two AllDifferent constraints in polynomial time. Our 
algorithm exploits a connection with matching on bipartite 
graphs. In particular, we consider simultaneous matchings. 
By generalizing Hall's theorem, we identify a necessary and 
sufficient condition for the existence of such a matching and 
show that the this problem is polynomial for convex graphs. 

Formal background 

Constraint programming. We use capitals for variables 
and lower case for values. Values range over 1 to d. We 
write D{X) for the domain of values for X, lh{X) {ub{X)) 
for the smallest (greatest) value in D{X). A global con- 
straint is one in which the number of variables ti is a pa- 
rameter For instance, AllDifferent([Xi, . . . , X„]) en- 
sures that Xi 7^ Xj for any i < j. Constraint solvers prune 
search by enforcing properties like domain consistency. A 
constraint is domain consistent (DC) iff when a variable is 
assigned any value in its domain, there are compatible values 
in the domains of all other variables. Such an assignment is 
a support. A constraint is bound consistent iBC) iff when a 
variable is assigned the minimum or maximum value in its 
domain, there are compatible values between the minimum 
and maximum domain value for all other variables. Such an 
assignment is a bound support. A constraint is bound disen- 
tailed iff no possible assignment is a bound support. 
Graph Tlieory. Solutions of AllDifferent correspond to 
matchings in a bipartite variable/value graph (Regin 19941. 

Definition 1. The graph G = {V,E) is bipartite ifV parti- 
tions into 2 classes, V — AU B and An B — such that 
every edge has ends in different classes. 

Definition 2. Let G = {AU B,E) be a bipartite graph. 
A matching that covers A is a set of pairwise non-adjacent 
edges MCE such that every vertex from A is incident to 



exactly one edge from M. 

We will consider simultaneous matchings on bipartite 
graphs(siM-BM) (Kutz et al. 200 8 ). 

Definition 3. An overlapping bipartite graph is a bipartite 
graph G = {AU B,E) and two sets S and T such that 
A = SUT, AnB = (!),andSr\T ^(d. 

Definition 4. Let {AU B, E) and S, T be an overlapping 
bipartite graph. A simultaneous matching is a set of edges 
MCE such that M D {S x B) and M D {T x B) are 
matchings that cover S and T, respectively. 

In the following, we use the convention that a set of ver- 
tices P is a subset of the partition A. We write N{P) for the 
neighborhood of P, = F n (S" \ T), = Pn{T\S) 
and P^^ = P n 5 n P. SIM-BM problems frequently oc- 
cur in real world applications like production scheduling and 
timetabling. We introduce here a simple exam timetabhng 
problem that will serve as a running example. 

Running example. We have 7 exams offered over 5 days 
and 2 students. The first student has to take the first 5 exams 
and the second student has to take the last 5 exams. Due to 
the availability of examiners, not every exam is offered each 
day. For example, the first exam cannot be on the last day of 
the week. Only one exam can be sat each day. This problem 
can be encoded as a SIM-BM problem. A represents the ex- 
ams and contains 7 vertices Xi to Xj. B represents the days 
and contains the vertices 1 to 5. S — [Xi,X2, X3, X4, X5] 
and T = [X3, X4, X5, Xg, X7]. We connect vertices be- 
tween A and B to encode the availability restrictions of the 
examiners. The adjacency matrix of the graph is as follows: 





1 2 3 4 5 


A'' = S\T 

X2 


* * * * 


X3 

A^"^ = SnT Xi 


* * * 

^ Hf^ :lf^ ^ -i^ 


A- = T\S f 
Xr 


* * * * 



Finding a solution for this SIM-BM problem is equivalent 
to solving the timetabling problem. □ 

Simultaneous Bipartite Matching 

We now consider how to find a simultaneous matching. Un- 
fortunately, this problem is NP-complete in general (Kutz] 
\et al. 20^8 ). Our contribution here is to identify a neces- 
sary and sufficient condition for the existence of a simultane- 
ous matching based on an extension of Hall's theorem ( |Hall| 
[1935). We use this to show that a simultaneous matching on 
a convex bipartite graph can be found in polynomial time. 

In the following, let G"^^ be the subgraph of the overlap- 
ping bipartite graph G that is induced by choosing an edge 
(u, u) to be in the simultaneous matching. If w e A^'^ then 
^(ti v) ~ ^ ~ ^i- If ^ £ (^nd symmetrically if 
u e V) then G"(„ = {V-{u},E\{{u',v)\u' ^S)). 
If M is a SIM-BM in G"(„ then AfU{(u, v)} is a SIM-BM 
in G. Since the edge [u, v) is implied throughout, we write 
G' = C?'(„^). In addition, we write N'{P) ^ Ng'{P). 



Extension of Hall's Theorem 

Hall's theorem provides a necessary and sufficient condition 
for the existence of a perfect matching in a bipartite graph. 

Theorem 1 (Hall Condition ( |Hall 1935| l). Let G = 

{AU B, E) such that A Ci B = 0^ There exists a perfect 
matching iff\N{P)\ > \P\for P C A 

Interestingly we only need a small adjustment for simul- 
taneous matching. 

Theorem 2 (Simultaneous Hall Condition (SIM-HC)). Let 

G ~ {AU B, E) and sets S, T be an overlapping bipar- 
tite graph. There exists a SIM-BM, iff\N{P)\ + \N{P^) n 
7V(P^)| > \P\forP C A. 

Proof We prove SIM-HC by induction on When |^| = 
1, the statement holds. Let |A| = fc > 1. 

If A'^ = or A'^ = then SIM-HC reduces to the condi- 
tion of Hall's theorem and the statement is true for that rea- 
son. Hence, we assume A^ ^ and A^ 7^ 0. We show that 
there is an edge (u, v) that can be chosen for a simultaneous 
matching and the graph G'^^^ will satisfy SIM-HC. Fol- 
lowing (Diestel 2006 ), page 37, we consider two cases. The 
first case when all subsets of A satisfy the strict SIM-HC, 
namely, |iV(P)| + \N[P^)r\N{P'^)\ > \P\ and the second 
case when we have an equality. 

Case 1. Suppose \N{P)\ + \N{P^) D N{P^)\ > \P\ 
for all sets P C A. As A^ 7^ we select any edge 
{u,v),u e A^ and construct the graph G"(„ (the case 

u £ A^ is symmetric). For any set P C A \ {u] we con- 
sider two cases: either v ^ N{P) or v £ N{P). In the 
first case, the neighborhood of P is the same in G and G", 
so the SIM-HC holds for P. In the case that v G N{P), 
then either z; is a shared neighbor of P^ and P^, which 
means that N'{P^) D N'{P^) = N{P^) D N{P^) - 1 
but N'{P) = N{P) by construction, or is a neighbor 
of P^ but not of P^. Therefore N'{P) > N{P) - 1. 
But N'{P^) n N'{P^) = iV(P^) n N{P^) by construc- 
tion. In either case, \N'{P)\ + |iV'(P^) n N'{P^)\ > 
|iV(P)| + |iV(P'5) n iV(P^)| - 1 > |P| for any set P in 
G". By the inductive hypothesis there exists a simultaneous 
matching in it. 

Case 2. Suppose that there exists a set P C A such that 

\N{P)\ + \N{pS)nN{P^)\ = |P|. LetQ^ {A'UB',E') 
such that A' ^ A\P, B' ^ B \ (A^(P^) n N{P^)) and 

E' = {{u,v) £ En{A' X B') \ 

(w e A^' w ^ N{P) \ NiP'^)) A 

(w e A^' =^ N{P) \ iV(P^)) A 

(wgA^T' =^ v(^N{P)) } 

There exists a simultaneous matching in G — Q by the 
inductive hypothesis. We claim that the SIM-HC holds also 
for Q. This implies that, by the inductive hypothesis, there 
exists a simultaneous matching in Q. Suppose there exists a 
set P' C A' that violates the SIM-HC in Q. 

We denote as N{P) the neighborhood of P in G and 
Nq{P') as the neighborhood of P' in Q. We know that 
the sets P' and P are disjoint. We observe that N{P U 



P') = iV(P) U N{P') = N{P) U Nq{P'), because 
{N{P')\Nq{P')) C Ar(P) by construction of Q. More- 
over |7V(P) U Nq{P')\ = \N{P)\ + \Nq{P')\ ~ \N{P) n 

Nq{P')\ = \N{P)\ + Wq{p')\ - mP) n (a^q(p^') u 

Nq{P^') U iVQ(P'^^'))|. By construction of Q, we have 
that Nq{P^'^') n iV(P) = 0, Nq{P^') n 7V(P) = 
A^qCP'^') n N{P'^) and 7Vq(P^') n A^(P) = ArQ(P^') n 
A^(P^). Hence, \N{P)\J Nq(P')\ = |Ar(p)| + _ 
where D = (iV(P'5) n Nq{P'^')) U (A^(P'^) n 
7Vq(P'5')). Similai-ly, 7V(P^UP^') 7V(P^) U iVQ(P^') 
and 7V(P^ U P^') = 7V(P^) U Nq{P'^'). Therefore, 

|iv(p^up^')n^(^^up^')| = |(7V(P^)u7VQ(P'5'))n 
(iV(P^) u ^q(p^'))| < I A^(P^) n 7V(P^)| + |iVQ(P^') n 



\D\. Finally, we have that 



|iV(PUP')| + |A^(P^ UP^ ) n A^(P'^ UP^ )| 



< 



\N{P)\ + \Nq{P')\-\D\ + 

\N{pS) n 7V(p^)| + |iVQ(p^') n iVQ(P^')| + = 
|A^(p)| + |iV(p^)niV(p^)| + 

\NQ{P')\ + \NQ{P^')nNQ{P'^'\ < 

\P\ + \P'\ = \PUP'\. 

Hence P U P' violates the SIM-HC in G, a contradiction. 
Therefore, there exists a simultaneous matching in Q. 

Let M be a simultaneous matching in Q. For any edge 
(u,v) G M, we construct the graph R ^ {G — Q){u,v) 
and show that Nr{P*) = Ng{P*) and iVj^(P^*) n 
NB.iP^*) = N{P^*) n N{P^*) for any P* C P. 

Let (u, v) be an edge in M. By construction of Q, we 
have that v ^ N{P^*) n N{P'^*). Hence, the construction 
of R leaves the size of Nii{P^*) f] A^i?(P^*) the same as 
in G. Moreover, u £ A^"^' ^ v ^ N{P*), u e A^' ^ v ^ 
N{P^*), ueA^'^vi NiP'^*). 

Consider the remaining options for the edge (u,v). If 
u e A^' (u e A'^' is similar) then v can be in iV(P^*), 
so V is a shared vertex with N{P*). There are two cases 
to consider: v G N{P^*) \ N{P^^*) and v e 7V(P^*) n 
N{P^^*). In the first case, the construction of R leaves 
the size of Nii{P^ ) the same as N{P'^ ) because u E 
A^' and can share vertices with P^*. In the latter case, 
NniP'^^*) = iV(P^^*) \ {v}, but NniP^*) = N{P^*). 
Hence, Nr{P*) = N{P*) in both cases. Therefore, the 
SIM-HC holds for any P* C P. □ 

Running example. In our running example A = 

[Xi, X2, X3, X4, X5, Xg, X7], A^ — [Xi,X2], A^'^ = 
X4, X^] and A"'" = [Xg, Xj] . It is easy to check that 
the simultaneous Hall condition holds for all subsets of the 
partition A. 

Note that Theorem [2] does not give a polynomial time 
method to decide if a simultaneous matching exists. Veri- 
fying that SIM-HC holds requires checking the exponential 
number of subsets of A. 



Removing edges 

To build a propagator, we consider how to detect edges that 
cannot appear in any simultaneous matching. 

Definition 5. Let G = {A\J B,E) and sets S, T be an 
overlapping bipartite graph. A set P, P Q A, is 

a simultaneous Hall set iff 

\N{P)\ + \N{P^)f^N{P^)\ = |P|. 

an almost simultaneous Hall set iff 

\N{P)\ + \N{PS) n iV(P^)| - |P| + L 

a loose set iff 

\N{P)\ + |7V(P^) n N{P^)\ > \P\ + 2. 

Tlieorem 3. G = {AU B, E) and sets S, T be an overlap- 
ping bipartite graph. Each edge {u, v), u £ A and v £ B 
can be extended to a matching that covers S and T iff 

1. for each set P: 

(a) \N{P)\ + \N{P^) n iV(P^)| > |P| 

2. for each simultaneous Hall set P: 

(a) ifu i P then v ^ (A^(P^) n N{P^)) 

(b) ifu£S\{TU P) then v ^ (iV(P) \ N{P'^)) 

(c) ifueT\{SU P) then v ^ {N{P) \ N{P^)) 

(d) ifu e{SnT)\P then v ^ N{P) 

3. for each almost simultaneous Hall set P: 

(a) ifu e{SnT)\P then v ^ iV(P'5) n N{P^) 

Proof. Soundness. The soundness of Rule[Ta|follows from 
Theorem|2] Let (u, v) be an edge that we want to extend to 
a matching. Suppose that {u, v) violates one of the rules for 
a sim-HallSet or an a-sim-HallSet P in G. We show 
that if {u, v) is selected to be in a matching, then P fails 
SiM-HC in G', 



Rule 2a 



If {u, v) violates Rule|2aJfor a sim-HallSet 
P then |5V1;P^) n N'{P^)\ = \N{W) f] N{P'^)\ - 1 and 
N'{P) = N{P), so the SiM-HC is violated for P in G' . 

Rule|2bj If {u,v) violates Rule pb| for a sim-HallSet 
P then \W\P)\ = \N{P)\ - 1 and\N'{P'^) D N'{P^)\ = 
|iV(P'5) n N{P^) \ so the SiM-HC is violated for P in G". 



|iV(P)| - 1 so the SIM-HC is violated 



Rule |2ct Symmetric to Rule|2b 

Rule|2dj If (u, v) violates Rule 2d for a sim-HallSet 
P then \W{P)\ 
for P in G". 

Rule 3a: If {u,v 
a-sim-HallSet P then 

|iv'(p^) n iv'(p^)| 
\N'{P)\ ■ 



Rule 



3a 



for an 

- 1 and 

- 1, 



violates 
\N'{P)\ = \N{PJ\ 
- = |iV(P^) n N{P^)\ - i, so 

\N'{P'^)r]N'{P^)\ = |P|- land the SIM-HC 
is violated for P in G". 

Completeness. Second, we show that Rules [2a 
cot7iplete 



3a\ are 

u,v) 



We will show that we can use any edge 
in a maching by showing that the graph G'^^ satisfies the 
SIM-HC, thus has a SIM-BM. 

Suppose there is a set P that violates the SIM-HC in G' 
but not in G so that 



and 



|A^'(P)| + |A^'(P^)niV'(P^)| < |P| (1) 



|iv(p)| + |7V(p^) niv(p^)| > |p| (2) 



Note that Af'(P) = N{P)\{v} md N'{P'^)nN'{P^) = 
N{P^)r\N{P'^)\{v}. Hence, \N'{P)\ > |iV(P)|-land 
|Af'(P^) n iV'(P^)| > |iV(P^) n N{P'^)\ - 1. 

There are three cases to consider for P in G, when P is 
a loose set, a sim-HallSet and an a-sim-HallSet in 
G. These cases are similar, so we consider only the most 
difficult case. Let P be an a-sim-HallSet in G. If 
u G then v (i N{P^) D 7V(P^) by Rule|3al Hence 



iV'(P'^) n N'iP' ) = 7V(P*) n N{P' ), 



so 



r'(p)i 



|iV'(P^')n7V'(P'^^ > \N(P)\ + \N{P'^)nN{P^)\-l> 
|P| and therefore (fTi and ^ cannot both be true. 

If u e {u e is symmetric) then v S N{P^) n 
7V(P^) or its complement. In the first case N'{P) = N{P), 
while in the second iV'(P^) n 7V'(P^) = iV(P'^) n 
7V(P''^). In both cases, |7V'(P)| + |7V'(Pf) n iVTP^)! > 
|A^(P)| + |iV(P^)n7V(P^)|-l > |P| so ([l) and (U cannot 
both be true. □ 

Running example. Consider again our running example. 
We show that Rules 2apa remove every edge that can 
not be extended to a matching. Consider the set P = 
{X2,X^,X4^,Xq}. This is a sim-HallSet as N{P) = 
{1,2,3}, N{P^) n 7V(P'^) = {2} and 4 = |P| 



2d 



\N{P)\ + \N{P^) n N[P'^)\ = 4. Hence, by Rule 

prune 1, 2, 3 from X^ and by RuleH^we prune 2 from Xi 
and Xj. Now consider the set P — {X2, X^, X^}. This is 



an A-SIM-HallSet. By Rule 3a we prune 2 from X4. The 
set P ~ {X2, X4, Xq} is also an A-SIM-HallSet and, 
by Rule [io] we prune 2 from X^. Next consider the set 
P = {X3,X4} which is a SIM-HallSet. By Rules [2fc 
and 2c we prune 1, 3 from Xi, X2, Xq and Xj. Now, {Xi\ 
is a SIM-HallSet and 4 is pruned from Xc, by Rule^^ Fi- 
nally, from the simultaneous Hall set {Xc,}, we prune 5 from 
Xi using Rule W^and we are now at the fixpoint. 





1 2 3 4 5 


= S\T 


Xi 




X2 






X3 




= S n T 


Xi 






X5 




= T\S 


Xe 




X-r 





The overlapping AllDifferent constraint 

We now uses these results to build a propagator. 
Definition 6. OverlappingAllDiff([X], S", T) where 

S C X,T C X, SUT = X holds iJALLDlFFERENT(5) 

and ALLDlFFERENT(r) hold simultaneously . 

Enforcing DC on the Overlap pingAllDiff constraint 
is A^P-hard (B essiere et al. 2007) 1. We consider instead en- 
forcing just BC. This relaxation is equivalent to the simul- 
taneous matching problem on a bipartite convex variable- 
value graph. Our main result is an algorithm that enforces 
BC on the OverlappingAllDiff constraint in 0{nd?) 
time. The algorithm is based on the decomposition of the 
OverlappingAllDiff constraint into a set of ai-ithmetic 
constraints derived from Rules [2b]-pal It is inspired by a de- 
composition of AllDifferent ( [Bessiere et al. 2009| |. As 



there, we introduce Boolean variables auu, bu to represent 
whether Xi takes a value in the interval [I, u] and the vari- 
ables C^ , C^^ and to represent bounds on the number 
of variables from S" \ T, T \ S" and S" n P that may take val- 
ues in the interval [l,u]. We introduce the following set of 
constraints for l<i<n, l<l<u<d and u — I < n: 



hi 


= 1 ^ X,<1 


(3) 




^1) = A b,u = 1) 


(4) 






(5) 












(6) 




ieS\T 








(7) 




ieT\S 




riST 




(8) 


riST 


+ Cl<u-l + l 


(9) 


riST 


+ Cl<u-l + l 


(10) 



We also introduce a dummy variable Cf^ = to simplify 
the following lemma and theorems. 

Lemma 1. Consider a sequence of values vi,V2, ■ ■ ■ ,Vk. 
Enforcing BC on ^ ensures u6(Cf„"^_]^) < 



<Vi+l 



Proof. For every i such that Vi < Vi+i, constraint (l8]l en- 
sures that ^&(Cf,^^,_i) < ub{C^J^_,) + ub{C^l^^_,). 
For every i such that Vi > u^+i constraint ([8]) ensures 
that u6(Cf„^^^_i) < «fe(Cf.T-i) - HCZ„..^i)- The 
left side of each inequality can be substituted into the right 
side of another inequality until one obtains ub{CfJ^_i) < 

Tlieorem 4. Enforcing BC on (|3]l-( |10| l detects bound dis- 
entailment 0/ OverlappingAllDiff in 0{nd^) time but 
does not enforce BC on OverlappingAllDiff. 

Proof. First we derive useful upper bounds for the vari- 
ables C^J . Consider a set P and an interval [a, b] such 
that N{P) = [a,b]. Let [ci,di] U ... U [ck.dk] be a 
set of intervals that tightly contain variables from P^ so 
that V^,[c„d,] e N{PS), [ei,/i] U ... U [e,„,/,„] be a 
set of intervals that tightly contain variables from P^ so 
that Vi, [ei,/i] e iV(P^), and /i U . . . U /p are inter- 
section intervals between intervals [ci,di] and [ei,/i], i.e., 

Vi,/, e iv(p^) n iv(p^). 

We first remove all intervals [ci^di] ([ei,/i]) that are 
completely inside an interval [ej,/j] {[cj^dj]). For any 
of these intervals [ci, di] ([e^, fi] is similar) there exists an 
intersection interval Ij such that Ij = [ci,di\. We de- 
note the set of removed intervals RI. The remaining in- 
tervals are {[ci, di], . . . , [cfe- , dfc']}, [d, /i] . . . , [e.^- , /„']} 
and Ii, . . . ,Ip'. For any interval [ci, di], ^ ensures that 
ub{Cl^.) <di-Ci + l- lb{C^.^J. Similarly for an inter- 
val [e-, /,], we have w6(CfJJ < - e, + 1 - lb{ClfJ. 



We sort the union of remaining intervals [ci,di] and 
[ci, fi] by their lower bounds and list them as semi- 
open intei-vals [31,32), [.93,54), [5fc'+m'-i>5fc'+m')- 
Using the sequence a, (51, 52, ■ • • , fffe'+m', 6 + l,a)^ 
where (.)^ indicates a repetition of x times the 
same sequence, Lemma [T] provides the inequality 

ubiC^ri) < ^&(C^fJ-i)+^("&(CfI-i) + ^KC,^J,,„,J + 

EtiubiCll) + Y.t,nb{Cl^f^) - IbiC^J)). Substi- 
tuting the inequalities that we already defined, we obtain 

ubiCf^,) < a-l + x{b~a + l + j:t, HCl,;) - 

YZ,Hclf;) + T.U\ii\-Hc!^))- 

For any removed interval Ij € RI we have — 

^KC^™(/^-)ma.(7,))) > Oor(|/,|-;6(C^,„(,^.)„^^(,^.))) > 

0. We reintegrate all removed intervals into the inequa- 



1 



x{b 



1 



tion to get u6(Cf'/_i) < 

Ell Hcl,:) - 1 Hcl,) + ELi \h\ - Hc^J))- 

Note that ELi l^.l = l^(^^) n N{P^)\, lb{C^,) ^ 
Ell HC^) IbiC^b) = Er=i 'K^^e^;/,)- Hence 
ub{C^D <a-l + x{b^a + l- IbiC^,) - lb{Cj,) 

+\NiP^)nNiP^)\~ibic^;f)) (*) 

Bound disentailment. Suppose, for the purpose of con- 
tradiction, that OverlappingAllDiff is bound disen- 
tailed and that constraints ((3|l-([T0l) are bound consistent. 
Then, there exists a set P, such that N{P) is an interval and 

|iV(P)| + |A^(P^)niV(P^)| < |P|. AsPfailsSlM-HC,it 
holds that lb{C^;f) + IbiC^J + /6(CJJ > |P| > |iV(P)| + 
|7V(P^) n A^(P^)| = 6 - a + 1 + |Af(P^) n iV(P^)| or 
lb{C!^) > b-a+2^lb{C!,)-lb{Cl,) + \N[P^)m{P^)\. 
Substituing the last inequality in (*) gives ub{CfJ_i) < 
a — 1 ~ X. Choosing a large enough value for x (say a) 
gives the contradiction ub{Ci^_i) < 0. 

Bound consistency. To show that this decompo- 
sition does not enforce BC, consider the conjunction 
of AllDifferent ([Xi,X2,X3]) and AllDifferent 
([X2,X3,X4]) with i?(Xi) - [2,3], D{X2) - [2,4], 
DiX^) = [1, 3], D{X4) = [1, 2]. Enforcing BC on 
does not remove the bound inconsistent value X2 — 2. 

Complexity. There are 0{nd) constraints (|3]l that can 
be invoked 0{d) times at most. There are 0{nd^) con- 
straints Q that can be invoked 0(1) times. There are 0{d^) 
constraints (|8]l that can be invoked 0{n) times. There are 
0{d^) constraints (|5]l- (|7]i that can be propagated in 0{n). 
The remaining constraints take 0{nd^) to propagate. The 
total time complexity is 0(ri(i^). □ 

It follows immediately that the simultaneous matching 
problem is polynomial on bipartite convex graphs. 
Theorem 5. A simultaneous matching can be found in poly- 
nomial time on an overlapping convex bipartite graph. 

Next, we present an algorithm to enforce BC. We show 
that constraints ([3]l-([T0]l together with the following two 
constraints enforce all but one of the rules from Theorem[3] 

c?: - c?; + i<i<u<d (11) 

lu — ^li 



Theorem 6. Constraints ((3|l-([T2|) enforce /?M/e^[2a]-[22/| 

Proof Sketch. Based on Lemma [T] similar to the proof of 
Theorem]?] we show that all intervals that contain variables 
from a SIM-HallSet P become saturated intervals, so that 
the lower bounds of the corresponding variables C'^^, C^ 
and C^ equal to their upper bounds. Hence, these values are 
pruned from domains of variables outside the set P. □ 

Theorem 7. Suppose constraints ^- \i2) together with 
2<l<k<u<d have reached 



r<ST _ r<ST 
'-'lu — '-'Ik 



their fixpoint. Rule 3a can now be enforced in 0{nd'^) time. 



Proof Sketch. Let P is an A-SIM-HallSet, A^(P) = [a, 6]. 
Similar to the proof of Theorem ]4] we can obtain that 
lb{C^^) + 1 > ub{C^^). Hence, we can identify inter- 
vals, that might contain an A-SIM-HallSet P. Next, we 
observe that if we add a dummy variable Z, D{Z) = [a, b] 
to the set P so that P' = P U {Z}, Z £ P-^^' then P' is a 
SIM-HallSet. This allows us to identify the set N{P^) n 
N{P'^) by simulating constraints (]3|)-(fT2]) inside the inter- 
val [a, b] taking into account the variable Z. There are 
0(d2) intervals. Finding N{P^) D iV(P^) takes 0(n + d) 
time inside an interval. Enforcing the rule takes 0{nd) time. 
Hence, the total time complexity is 0{nd'^). □ 

From Theorems ]6] and]?] it follows that 

Theorem 8. BC on OverlappingAllDiff can be en- 
forced in 0{nd^) time. 

Running example. We demonstrate the action of 
constraints (]3|)-(]T2]). The interval [1,4] contains a 

SIM-HALLSET P ^ {X2,X3,X4,X(i}. 

lb{C^2) >±and lb{C^r^)_^ 1 ^ ^ 



2d 



and ub{C^^ ) < 1 ^ (|8j 



ub(Gft^ 



Rule 

(]T0| ^ ub{Cff) <l 
) < 2. The interval 
ub{Cf^). Hence, 



l<l<u< d 



(12) 



[1,3] is saturated, as lb{Cf^) 

by (]3])-(|5l [1,3] is removed from DiX^). "Rule s \2b\2c. 
As lb(d^) = 2 ^ ^ ^ "^(Cfs) < 1 ^ O % 
ub{Ci2) ^ 1- The interval [1,2] is saturated, as 
lb{C^2) = ub{C^^). Hence, by [1, 2] is removed 

from D{Xi). Similarly, [2,3] is removed from D{Xi). 
Rule^a^ This is satisfied as 2 is removed from all variables 
outside P. 

Exponential separation 

We now give a pathological problem on which our new prop- 
agator does exponentially less work than existing methods. 

Theorem 9. There exists a class of problems such that en- 
forcing BC on OverlappingAllDiff immediately de- 
tects unsatisfiability while a search method that enforces DC 
on the decomposition into ALLDIFFERENT constraints ex- 
plores an exponential search tree regardless of branching. 

Proof. The instance Z„ is defined as follows I„ = 

ALLDIFFERENT([X U Y]) a ALLDlFFERENT([y U Z]), 
D{Xi) = [1, 2n - 1], i = 1, . . . , n, D{Y,) = [1, An - 1], 
i = 1, . . . , 2n and D{Z,) = [2n, 4n - 1], z = 1, . . . , n. 

OverlappingAllDiff. Consider the interval [l,4n - 
1]. |P| = An, \N{P)\ = 4n-l and |A^(P^)n A^(P^)| = 0. 
By Theorem ]2j we detect unsatisfiability. 



Decomposition. Consider any AllDifferent con- 
straint. A subset of n or fewer variables has at least 2n — 1 
values in their domains and a subset of n + 1 to 3n variables 
has 471—1 values in their domains. Thus, to obtain a Hall set 
and prune, we must instantiate at least n — 1 variables. □ 

Experimental results 

To evaluate the performance of our decomposition we car- 
ried out an experiment on random problems. We used 
Hog 6.2 on an Intel Xeon 4 CPU, 2.0 GHz, 4GB RAM. 
We compare the performance of the DC, BC (L opez-j 



Ortiz et al. 2003) propagators and our decomposition 
into constraints ([3]i-([T2]i for the OverlappingAllDiff 
constraint (oBC). We use randomly generated problems 
with three global constraints: AllDifferent(X U W), 
AllDifferent(FUT/F) and AllDifferent(ZUVF), and 
a linear number of binary ordering relations between vari- 
ables in X, Y and Z. We use a random variable ordering 
and run each instance with 50 different seeds. As Table [T] 
shows, our decomposition reduces the search space signifi- 
cantly, is much faster and solves more instances overall. 



Table 1: Random problems, n is the size of X, Y and Z, 
o is the size of W, d is the size of variable domains. Num- 
ber of instances solved in 300 sec out of 50 runs / average 
backtracks/average time to solve. 



n,d,o 


BC 


DC 


OBC 




#s / #bt / t 


#s / #bt / t 


#s / #bt / t 


4, 15, 10 


14/2429411 /61.8 


41 /1491341 /52.1 


42/ 17240/ 32.5 


4, 16, 11 


6/5531047/ 153.7 


22/1745160/67.9 


31/ 8421/ 19.5 


4, 17, 12 


1/17/0 


6 /2590427 / 100.9 


24/ 8185/ 21.5 


5, 16, 10 


1 1 /3052298 / 82.0 


37/1434903/58.2 


42/ 20482/ 48.5 


5, 17, 11 


2/3309113/94.5 


19/2593819/ 114.6 


26/ 4374/ 15.8 


5, 18, 12 


1/17/0 


4/2666556/ 133.1 


22/ 3132/ 12.2 


6, 17, 10 


11 /2845367/79.1 


31 /143I671 /66.3 


40/ 6796/ 21.9 


6, 18, 11 


4/ 199357/6.6 


16/1498128/80.2 


31 / 4494 / 17.5 


6, 19, 12 


4/3183496/ 110.0 


5 /1035126/66.2 


27/ 3302/ 15.5 


TOTALS 








sol/total 


54 /450 


181/450 


285 /450 


avg time for sol 


78.072 


70.551 


24.689 


avg bt for sol 


2818926 


1666568 


9561 



Conclusions 

We have generalized Hall's theorem to simultaneous match- 
ings in a bipartite graph. This generalization suggests a 
polynomial time algorithm to find a simultaneous match- 
ing in a convex bipartite graph. We applied this to a prob- 
lem in constraint programming of propagating conjunctions 
of AllDifferent constraints. Initial experimental results 
suggest that reasoning about such conjunctions can signifi- 
cantly reduce the size of the explored search space. There 
are several avenues for future research. For example, the al- 
gorithmic techniques proposed in ( ,Puget 1998j ) and ( jLopez- 
Ortiz et al. 2003 1 may be generalizable to simultaneous bi- 
partite matchings, giving more efficient propagators. Fur- 
ther, matchings are used to propagate other constraints such 



as N Value ( Bessiere et al. 2006) l. It may be possible to ap- 
ply similar insights to develop propagators for conjunctions 
of other global constraints, or to improve existing propaga- 
tors for global constraints that decompose into overlapping 
consti-aints like SEQUENCE ( Brand et al. 2005) . Finally, 
we may be able to develop polynomial time propagators for 
otherwise intractable cases if certain parameters are fixed 
( ,Bessiere et al. 20 08 ) . 
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